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Abstract. We consider the reduced density matrix of a large block of consec- 
utive spins in the ground states of the XY spin chain on an infinite lattice. We 
derive the spectrum of the density matrix using the expression of the Renyi 
entropy in terms of modular functions. The eigenvalues A n form an exact geo- 
metric sequence. For example, for strong magnetic field A n = C exp (— nron), 
here to > and C > depend on the anisotropy and the magnetic field. Differ- 
ent eigenvalues are degenerated differently. The largest eigenvalue is unique, 
but the degeneracy g n increases sub-exponentially as eigenvalues diminish: 
g n ~ exp (ir^/ n/3). For weak magnetic field expressions are similar. 



1. Introduction 

Entanglement is a peculiar feature of a quantum system, which distinguishes 
it from a classical one. While the notion of entanglement has been introduced 
since the dawn of quantum mechanics, only recently physicists have employed a 
quantitative approach, mostly inspired by the progresses in quantum information 
and Bethe ansatz. 

The most studied quantity is the Von Neumann entropy, which is the quantum 
analog of the Shannon entropy and measures the amount of (quantum) information 
stored in a system A described by a density matrix p A : 

(1.1) S(p A ) = -Tr(p A lnp A ) . 

We shall consider the simplest case of a pure system U described by a wave function 
The system U is composed by the union of two disjoint subsystems A and B, 
meaning U = A U B. The reduced density matrix of subsystem A is obtained by 
tracing out the B degrees of freedom p A = Trs|^')(^ r | and the Von Neumann en- 
tropy of a subsystem is a measure of the entanglement between the two subsystems 

m 
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Characterizing the entanglement with a single number is definitely appealing, 
but it can hardly capture its complexity. For this reason, other entanglement 
measures have been introduced, such as, the Renyi entropy [2l [3l [4l [5l [6] : 

(1.2) S R ( PA ,a) = -^—lnTr(pl), l>a>0. 

1 — a 

Note that in the limit a — > 1, the Renyi entropy recovers the Von Neumann entropy: 

lim S R (p Al a) = S(pa)- 

We calculated the Renyi entropy of a large block of spins for the XY spin chain 
in [111 [8] and studied its analytical continuation into the complex plane of the 
parameter a. 

Another quantity that recently has attracted a lot of interest is the spectrum 
of the reduced density matrix, which is now referred to in the literature as the 
entanglement spectrum, after [7]. The knowledge of the entanglement spectrum 
fixes the density matrix up to unitarity transformation, we can say that it fixes the 
state of the block. It is also clear that the Renyi entropy and the entanglement 
spectrum are also very closely related. In fact, if we know ^-function of pa- 

oo 

(1-3) Q PA ^)^tvp a A ^Y.9nK- 

n=0 

at all a then we can find the eigenvalues A„, which have the meaning of probabilities 
(0 < A„ < 1) and their multiplicities g n . From the definition of the Renyi entropy 
(|1.2[) we see that 

(1.4) Cpa(") =e (1 ^ a)5n(pA ' a) • 

The aim of this paper is to use this relation to calculate the spectrum of pa for 
the anisotropic XY model, for which analytical expressions for the Renyi entropy 
are known explicitly [8] . The XY model is one of the simplest integrable mod- 
els (since quasiparticle excitations are essentially free fermions), while still having 
an interesting and non-trivial phase diagram. Von Neumann entropy of the XY 
model was first determined rigorously using Toeplitz determinant representation 
and Riemann-Hilbert techniques in [9] and later using thermodynamic arguments 
in [10j . and it was studied in |12j . The method of [9] was generalized in [13] for 
more general quantum spins. The evaluation of the Renyi Entropy for the XY 
model, which was done in [8] , uses again the Fisher-Hartwig formulae and the 
Riemann-Hilbert approach. Although this is the first time that these exact results 
are used to access the full spectrum of p A , general behaviors have already been 
known, because of the underlying free fcrmionic structure. 

In fact, in [7], in analogy with the formulae of standard statistical mechanics, 
it has been proposed to represent the reduced density matrix as pa — g~ Ua jZ. 
For AKLT models of interacting spins this formula was rigorously proved and the 
entanglement spectrum was calculated for several of these systems (see [141 115j). 
For essentially non-interacting models, such as the XY spin chain, the entanglement 
spectrum is known to be equidistant and the degeneracy of the eigenvalues is given 
by the number of ways a given energy can be realized by different excitations. This 
is essentially a partitioning problem and it has been addressed already several years 
ago in an effort toward the optimization of Density Matrix Renormalization Group 
approaches |16j . Our exact approach will agree with the results of |16j . 
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Figure 1 . Phase diagram of the anisotropic XY model in a con- 
stant magnetic field (only 7 > and h > shown). The three 
cases 2, 1a, 1b, considered in this paper, are clearly marked. The 
critical phases (7 = 0, h < 2 and h = 2) are drawn in bold lines 
(red, online). The boundary between cases 1a and 1b, where the 
ground state is given by two degenerate product states, is shown 
as a dotted line (blue, online). The Ising case (7 = 1) is also 
indicated, as a dashed line. 



It should also be mentioned that for some spin models the reduced density 
matrix was evaluated explicitly. 



2. Quantum entropies for the XY model 

The anisotropic XY spin chain is defined by the Hamiltonian 

00 

(2.1) H=- J2 [(1 + 7K^ +1 + (1 - 7VM+i + K] > 

J = -O0 

where 7 is the anisotropy parameter, erj, en and cr| are the Pauli matrices and h is 
the magnetic field. The Hamiltonian is clearly symmetric under the transformations 
7— > — 7 or ft - > — h, therefore we can consider just the quadrant 7 > and h > 
without loss of generality. The system is gapped in the bulk of the phase diagram 
and has two phase transitions where the spectrum becomes critical: at 7 = 0, \h\ < 2 
one has the XX model (universality of free fermions on a lattice) and \h\ = 2 is 
the critical magnetic field of the Ising phase transition, see Figure [TJ 
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The model was solved in |19|, I20|, 121), 122] . It is known that its correlation 
functions can be calculated using methods of Toeplitz determinants and Ricmann- 
Hilbert problems. These techniques were applied in [9l 1111 1121 [8] to the study of 
the quantum entropies. 

The density matrix of the unique ground state \GS) of the model is given by 
Pab — \GS)(GS\. The reduced density matrix of a subsystem A is pa = Ttb(pab)- 
We take the subsystem A to be a block of n consecutive spins (system B is the state 
of the rest of the chain) and consider the double scaling limit l<<n<<7V = oo, 
where N is the total number of sites in the chain, which we take to be infinite. 
Matrix elements of pa are correlation functions, see formulae (17) and (18) of |11) . 

The Renyi entropy (|1.2|) converges to the von Neumann entropy (|1.1|) for a — > 1, 
therefore we can concentrate just on the former quantity. Its analytical expressions 
were derived in |5j and they can be written as 
(2.2) 

1 a 1 1 , / 0f(O|mr o ) \ 1 



S R ( PA ,a) = - ln(fe k') + - In - ^ 1 u . + -In 2 , 

61 — a 31 — a \y 2 (0|air ) 04(Q\aiT O ) / 3 

for h > 2 and 
(2.3) 

S R (p A ,a) = -- In - +-- In 21 . +-ln2, 

6 1 — a \K / 3 1 — a \6s(0\atTo) 6i(u\a%To) J 3 



for h < 2; where 

(2.4) k' = V^~P, 

I(k) 

I(k) is the complete elliptic integral of the first kind, 

dx 



(2.5) I(k) = 



l ^/{l-x 2 )(l-k 2 x 2 ) 
and 

(2.6) 9j(z\T):=9j(z,q), q = <? iT j= 1,2,3,4, 

are the elliptic theta functions defined by the following Fourier series (|g| < 1) 

-1 \ 2 



(2.7) O x {z,q) = i (-l) m 9 (22Fi) G M m -h) 

771= — OO 
OO 

(2.8) e 2 (z,q) = J2 q(^ 2 e 2i < m -^ , 



771= — OO 
OO 



(2.9) e 3 (z,q) = 53 q"e 2izm 



TIT, — — OO 



(2.10) 4 (^<Z) = J2 (- 1 ) m « 



m „m g2izm 
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The elliptic parameter k = k("j,h) is defined in the different regions of the phase 
diagram as 



(2.11) k 



vW2) 2 +7 2 -l/7 



Case la: 4(1 - 7 2 ) < h 2 < 4; 



v /(l -/i 2 /4-7 2 )/(l -fe 2 /4) , Case lb: h 2 < 4(1- 7 2 ); 
7 / y/Qi/2) 2 + 7 2 - 1 , Case 2 : /i > 2. 



Alternatively, we can write the Renyi in terms of the A - modular function (see 
[§]) as 
(2-12) 



Sr = 



' " In(fcfc')- 1 ' 



In 



li 1 a 12 1 - a 

6 1 - a Vfc 2 / 12 1 - a 



A(air ) (1 - A(mr )) 
1 — A (ad to) 



A 2 (air ) 



iln 2 

3 



- In 2 ft > 2 
3 

h< 2 



The modular function is defined as (9r > 0) 



(2.13) 



A(r) = 



0f(O|r) 



T W(r), and 1 - A(r) = gffi = , 



9l(0\r) 



We note that the basic modular parameter k = k{"y, h) defined in (12.111) coincide 
with the value of the function k(r) at r = ir . 

A third representation of the Renyi entropy in terms of q— series will be useful 
to determine the multiplicities of the reduced density matrix eigenvalues 
(2.14) 



S R (fl A ,a) 



1 a 
12 1 - a 
1 a 



In 



6 1 -a 



fc 2 fc' 2 
Wq 
, /16#' 



1 — a 
2 

1 -a 



In 



n[ 

m=0 

OC 

InJ] [1 



1 



g (2ra+1)a ft > 2 
? 2ma ]+ln2 ft< 2 



where 
(2.15) 



--nl(k')/l{k) 



3. Spectrum of pa 

Using the expressions for the Renyi entropy we just listed, we now want to de- 
termine the eigenvalues A„ (0 < A„ < 1) of the operator pa and their multiplicities 
g n , through its momentum function (ll.3[) using (|1.4[) . Using (|2.14l) we have 



(3.1) 

where 
(3.2) 



no. j-i iti kJsL 



2e 



6 6 4fc. 



n( 1 + ^ m+1 ) 2 h>2 

■0 

00 

IK 1 



m=0 



g 2m ) 2 ft<2 



m— 1 



g a = e- Q7rro = q Q 



To use these expression, we will need some results on q-series and elementary notions 
of the theory of partitions. 
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Let us concentrate first on the case h > 2. Classical arguments of the theory 
of partitions (see e.g. [23]) tell us that 

oo oo 

(3.3) n( i + « 2n+i )=E^ i) H'z n , 

n=0 n=0 

where p@ (0) = 1 and p^(n), for n > 1, denotes the number of partitions of n into 
distinct positive odd integers, i.e. 

(3.4) p£\n) = #{K, . . . ,m k ) : m 3 ■ = 2rj + 1, mi > m 2 > . . . > m k , 

n = mi + m-x + • • • + m k } ■ 
Hence (|3.1[) for h > 2 becomes 

(3-5) ( PA (a) = e a ^ + -° ln ^)J2 a n<l2, 

n=0 

where, 

n 

(3.6) a = 1 , a n = (QPo( n ~ 
Since 

(3.7) <£ = (e- WT °") Q , 
we conclude that 

oo 

(3.8) C P .(«)=E a » A - A„=e-^«+^+^n^. 

n=0 

Comparing the last equation with (|1.3p we arrive at the following theorem: 

Theorem 3.1. Let the magnetic field h > 2. Then, the eigenvalues of the reduced 
density matrix pa are given by 

1 i n k k' J(fc') r„ 11 

(3.9) A n = e sln ~" 7r wi n -raJ, n = 0,1,2,... 

and the corresponding multiplicities g n — a n are defined by (3. 



The case h < 2 is treated in a very similar way. Instead of Q3.3P we use another 
combinatorial identity, 



(3.io) nu+^HEffiW". 



»=1 n=0 

where pV/"* 1 (0) = 1 and pf^ (n), for n > 1, denotes the number of partitions of n into 
distinct positive integers, i.e. 

(3.11) P { ^{n) = #{(mi, . . . ,mfc) : mi > m 2 > ■ ■ ■ > m k > 0, 

n = mi+m 2 + ... + m k }. 
It is worth noticing that (see e.g. |23j ) 

(3.12) p$(n)= Po (n), 

where pe>(n) denotes the partitions of n into positive odd integers: 

(3.13) po(n) = #{ (mi, ... , m fc ) : mj = 2r,- + 1, n = mi + m 2 + . . . + m fe }. 
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The analog of equation (|3.5[) , with the help of (|3.10|) , now reads as 

(3.14) C PA (a) = 2e Q (-^+^)£6 ng - 

oc 

= 2j2b n \«, A„ = e- 2 — -^H^, 

n—0 

where 

n 

(3.15) 6 = 1 , &„ = 5>#(0P#(n - • 

Finally, comparing (|3.14l) with equation (|1 .3[) we arrive at the analog of Theo- 
rem [34] for the case h < 2: 



Theorem 3.2. Lei the magnetic field h < 2. Then, the eigenvalues of the reduced 
density matrix pa are given by the equation, 

(3.16) A„ = e s ln ^- 27r w["+A], n = 0,1,2,... 

and £/ie corresponding multiplicities g n — 2b n where the integers b n are determined 

by cm. 

4. Asymptotics of g n 

Consider first the case h > 2. Following the usual methodology, we introduce 
the generating function 



(4-1) /(*):= X>„; 



n=0 

This function is holomorphic in the unit disc. Indeed, we have from (I3.5[) that, 

(4.2) f(z)=e- a ^ in ^ + ^)( PA (a), a = lnz. 

7TT 

Statement of holomorphicity then follows from the first equation in (|3.ip . 

The function f(z) has a singularity at z = 1. In order to sec this, we deduce 
from (|1.4|) the representation for f(z) in terms of the entropy Sr{pa, ot) 

(4.3) /(*) = e ^)s R(PA , a) - a (^^ + ^)^ a = _^ lnz _ 

7TT 

In [5], using the explicit formulae (|2.2I) and (|2.12p . and the modular properties of 
the A - function, 

(4.4) =1-A(r), A(r + 1)- A(rl 
it was obtained that 



A(r)-1' 



(4 ' 5) ^' a) = ^i^)12 7M + l-^6 ln — + °( e "J 



a(l-a) 12 J(fc') 



7T 7T 

a — » 0, < arga < — . 

2 6 2 



Hence, 

(4.6) f[z) = e _T ^ + ^ ln " +0 ( c " 2/ln "), z->l, \z\<l 



F. FRANCHINI, A. R. ITS, V. E. KOREPIN, AND L. A. TAKHTAJAN 

Rs(log(f(p„»'"))-Log»(p„» r ")) 



Lo g [f(,»|-nLogp 





Figure 2. Case h > 2: plots of the logarithm of the generating 
function in the Cauchy integral (|4.7j) as a function of the radius p 
(left panel) and of the angular phase 9 at the saddle point radius p n 
given by (|4.10j) (right panel; only real part shown). The plots show 
the comparison between the exact expression (|4.3I) (continuous line, 
blue on-line) and its asymptotic approximation (14.61) (dashed line, 
red on-line) at n — 40. 



(4.7) 



The coefficients g n are given by the Cauchy formula, 

2tti y w=1 _ e 



From this formula, the large n asymptotics of g n can be rigorously obtained by 
the Hardy- Ramanuj an- Rademacher circle method (see e.g. 125 ) using the modular 
properties (|4.4[) of the A-function. According to the circle method, the leading 
contribution to integral (14.7[) comes from the neighborhood^ of the point 1 — e. 
This fact can be also demonstrated by plotting the function g(z) = In f(z) — n In z, 
see Figure [2] Therefore, we can replace the explicit formula (|4.7j) by the estimate, 



1 



,n+l 2 7ri 



-GO) 



Sn 2tti J :"" 1 I'm ./^ 

where £ = {|z| = p„, |argz| < <5 < — lnp„}, and 

1 



(4.9) 



G(z) 



n In z 

12 



dz 



121nz 

We determine p n as the stationary point of G(z), i.e.: 
dG{z) 

z=p 



(4.10) 



dz 







n — > oo, 



The implementation of the circle method in its full power would yield the Hardy-Ramanujan- 
Rademacher type expansions for the multiplicities g n (cf. 25 where the classical case of pjv(n) 
is considered). In this paper, we are only concerned with the leading behavior of g n , and to this 
end we only need the localization of the integral near the point z = 1 — e. The rigorous proof of 
this property of integral 14.71 is not at all trivial. Indeed, it needs again the modular properties 
of the function /(z) . 

It should be also mentioned that there are more general techniques of the asymptotic analysis 
of the partitions, such as Meinardus theorem (see e.g. |23| ; see also 24). These techniques do 
not exploit the modular properties of the corresponding generating functions, however, unlike the 
circle method, they only provide the leading terms of the asymptotics. The Meinardus theorem, 
as it is stated in |23| , is not directly applicable to generating function 14.111 . 
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Switching then to polar coordinate z = e e ™+ ie , we can rewrite G(z) in the form, 



(4.11) 



i 



12 12d 



1 + 



z€ C. 



It can be shown, using again the circle method, that there exists a positive kq 
such that the following choice of the parameter 5 in the definition of the arc L is 
consistence with estimate 



(4.12) S = el+ 2K , < k < Ko- 

Using this specification of S we re-write estimate (|4.11l) as 

^2 



1 + O e 



(4.13) 



1 + (n K )} , < k < n , n -> oo, ze£. 



V3V'^ 12 ' 12e3' 
In its turn, this estimate yields the following representation for G(z) on £ 



(4.14) 



G(z) = -^=Jn ' 



12 12e; 



where is a function holomorphic in the neighborhood of the interval [—8, 5] and 
satisfying the estimates, 



(4.15) 



^)=0(l + O(n- K )), |^l + 0(n- K ), 



0<k<k o , n — > oo, e [—6, 5}. 
From (|4. 14[) we have that 



(4.16) 



2vri 



-G(z) 



• dz = e 



2tt 



-t 2 («) 



d0. 



—5 



At the same time, equations (|4.15j) allow us to use t = t(9) as a new integration 
variable and transform (|4.16[) into the asymptotic formula, 



27ri 



-GO) 



■ dz = e Va 



2tt 



t(S) 



t{-5) 



t 2 d9 
~dt 



dt 



: 



V n -T2 



2n 



— dt + O [n 2 e " 
dt V 



(4.17) 



e 75 



2tt 



2 2 

^* dtfl + O + O (n- 1 ) 



n — > oo, 
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Assume now that 



(4.18) 0<K<min|i,K O 
Then, the integral in the right hand side of (I4.17|) can be estimated as follows, 
— / e ™% dt= — e ^ dt / e ^ di 

If 00 x 2 ,2 , 1 /-OO 2 1 /-OO „2 2 

= J_/ e -i^* d ^_ e 3/2i / e -^ dt= l / e "^* di + Ofn"" 



(4.19) 



J-a/H^ + O (n~~) = 2- 3 / 2 3- 1 / 4 fn - ^) ^ + O (n^) 



Estimates (|4.8|l . (|4.17p . and (|4.19|l yield the following asymptotic formula for the 
multiplicities g n of the eigenvalues of the reduced density matrix for h > 2. 

(4.20) g n ~ 2- 3 / 2 3- 1 / 4 n- 3 / 4 e 7r v / t , „ oo . 

Turning now to the h < 2 case, using f|3 . 14[) and remembering that g n = 2b n , 
we have 

(4.21) Q PA {a) = e°(* ln ^"^) / (e" 2 ^") , 

with /(z) defined as in (|4.ip . As before, using (|1 .4[) we can express /(z) in terms 
of the Renyi Entropy Sr(pa, a): 

(4.22) / (^ = e a-«)SR(w,a)-a(iln^-^) ) q= L^. 

27TTO 

Again, we are interested in the neighbors of the point z ~ 1 (see Figure [3]), 
where we can use the asymptotics derived in [8]: 



(4 - 23) ^' a) = a<1^12 7M + T^6 ln 4P + °( e 

7T 7T 

a^O, --<arga<-. 

Using ((47231) in (t4T22|) we have (cf. (l4~6]) ) 

(4.24) /(z) = e"^" Alnz+ °(° 2 " 2/lni ), |z| < 1 . 

The calculation proceeds exactly as before, where in (|4.8I) we now have 



(4.25) G( ^__ + n+ )ln , 



and the saddle point z = p n is: 
(4.26) ^ 



=>■ Pn = e v^STTT = e 



^It is worth noticing, that under condition 14. 181 the term O (n 1 ) in 114. 1711 becomes in fact 
O (e-CK™' 9 ), > , that is O (n" 00 ). 
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Re (log(f( p„,'> )) - Log n ( p„«" )) 



L.,»[f(,»|-nLogp 





Figure 3. Case h < 2: plots of the logarithm of the generating 
function in the Cauchy integral (|4.7j) as a function of the radius p 
(left panel) and of the angular phase 9 at the saddle point radius p n 
given by (|4.26j) (right panel; only real part shown). The plots show 
the comparison between the exact expression (I4.22[) (continuous 
line, blue on-line) and its asymptotic approximation (|4.24l) (dashed 
line, red on-line) at n = 40. 



Instead of ()4.17|) we now have, 
(4.27) 

J_ r e- G « 

2tti 



dz = C 7r V / f(™+T2 



1 

2^ 



Wl + O (n - *)) + 0(n" 1 ) 



with S defined by exactly the same equation f|4. 12|) . Choosing k as in (|4.18p . we 
can again approximate the integral in the right hand side of (I4.27|) by a complete 
Gaussian integral (cf. 14.191) : 



(4.28) 



1 

2^ 



dt 



1 

2^ 



1 /fir 3 

±J™H + („-<*>) 
2n V 7r v ' 



2-5/43-1/4 f n + _L 



-3/4 



0{n 



Estimates (|4~8|) . (|4.27[) . and (|4.28[) yield the following asymptotic formulcQ for the 
multiplicities g n of the eigenvalues of the reduced density matrix for h < 2. 



(4.29) 



9n 



2 -5/4 3 -l/4 n -3/ Vv /l 



n — ^ 00 



5. Critical lines 

As we mentioned in the introduction, there are two critical lines in the phase 
diagram of the XY model, where the gap closes: at the critical magnetic field 
h — 2 and at the isotropic line 7 = 00 and \h\ < 2. In critical cases, the £ function 
behaves as 



(5.1) 



4 It is worth noticing that for the case h < 2 the generating function f(z) can be easily 
transformed to the one satisfying the conditions of the Menardus theorem and hence asymptotics 
1 14. 291 1 can be also obtained by using the Menardus theorem. 

^Also known as the XX spin chain. 
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where T a is a non-universal constant and £ is the relevant length-scale in the con- 
sidered regime. This was first discovered for the isotropic case [11] and confirmed 
by conformal field theory |26j . The straightforward application of the asymptotic 
formulae found in [8J, shows agreement with (|5.1I) close to the critical lines, with 
the length-scale set by the inverse energy gap A, i.e. 

(5.2) f" 1 =A = \h-2\ , c=i, h->2, 

(5.3) r 1 = A = N, c = l, -f->0,\h\ <2. 

Our Theorems 13.11 and 13.21 for the eigenvalues distributions in the different 
regimes remain valid arbitrarily close to the critical lines: while the eigenvalues X n 
tend to collapse and vanish with the energy gap, namely 

(5.4) A n ~ At e^™ , A -> , 

their multiplicities do not depend on the parameters 7 and h of the model. This 
agrees with [28] . 

6. Conclusions 

We have calculated the spectrum of reduced density matrix of a large block of 
spins in the ground state of XY spin chain [entanglement spectrum [7] ] , using our 
results on Renyi entropy [§]. 

We have confirmed the expectation that, being the model essentially non- 
interacting, the eigenvalues are equidistant and their multiplicities have simple 
interpretation in terms of combinatorics and different partitions of integers, see 
Theorem 13.11 and 13.21 The exact formulae for the eigenvalues have been given in 
terms of the parameters of the model. 

The asymptotic behavior of the multiplicities has been calculated, using the 
modular properties of the Renyi entropy. The leading terms of the asymptotics are 
given in equations (|4.20l) and (|4.29[) for strong and weak magnetic field, respectively. 
Our results agree with the estimates in |16j . However, this is not the log- normal 
behavior quoted in |18j for the XY model and also taken from |16j . In fact, the 
log-normal result was achieved by combining (and smearing) the degeneracy with 
the eigenvalue behavior to give an estimate of the behavior of an effective eigenvalue 
in a non-integrable system. This estimate is important to implement an efficient 
DMRG calculation for generic systems. 
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